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ABSTRACT 

A critical- speed analysis of undamped rigid rotors on flexible foundations reveals 
four conical- mode solutions comprising forward and backward precession. There is also 
a two- solution bouncing- mode set dependent only on foundation- to -rotor mass and spring- 
constant ratios. Although analysis of a firm-foundation model omits two of the four 
conical- mode solutions, good agreement exists between the two corresponding solutions. 
The flexible-foundation analysis can be approximated by a simple firm-foundation result 
when the foundation- to- rotor spring- constant ratio is zero or infinite. Only three calcu- 
lations are needed in this range to cover all intermediate values. 
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CRITICAL-SPEED ANALYSIS OF RIGID ROTORS ON FLEXIBLE FOUNDATIONS* 

by Richard H. Cavicchi 
Lewis Research Center 

SUMMARY 

A theoretical analysis of rigid rotors in undamped bearings on flexible foundations 
develops a general frequency equation for both forward and backward precession. One 
set of two unique solutions comprises the bouncing- mode natural frequency, which is a 
function of foundation- to-rotor mass and spring- constant ratios, but independent of rotor 
speed and moments of inertia. Another set, the conical mode, is a function of rotor 
speed and moments of inertia and contains four unique solutions. Example maps are 
given for a wide range of rotor speed and of mass, moment- of-inertia ratios, and 
spring- constant ratios as a guide in preliminary design. Besides locating major critical 
speeds, the maps locate nonsynchronous critical speeds that may result from bearing de- 
fects. 

The four coni cal- mode critical- speed solutions appear on the maps as a high- and 
low-frequency set of two-branch curves*. Although the firm-foundation analysis omits 
two of these critical speeds, its two existing solutions agree well with flexible-foundation 
results. Furthermore, the two omitted critical speeds can be approximated by the two 
firm-foundation bouncing- mode solutions. The flexible-foundation analysis can be ap- 
proximated by the two- solution firm- foundation result for the zero and infinite extremi- 
ties of the foundation- to- rotor spring- constant ratio. Therefore, only three flexible- 
foundation calculations are needed within this range to cover all intermediate values. 

INTRODUCTION 

Stationary turbo machinery is customarily mounted on rigid foundations. Therefore, 
many theoretical studies of rotor vibration assume firm foundations, as in reference 1. 
An important extension is to regard the foundation itself as flexibly supported. Such con- 
figurations are more representative of aircraft and spacecraft installations. Analysis of 
models for this application is complicated by the addition of at least two variables to 
those of firm-foundation models. 

*The original version of this report has been corrected to use gyroscopic terms 
based on the rotor absolute angle instead of on the relative angle between the rotor and 
foundation. 



Crook and Grantham report in reference 2 on a method of predicting synchronous 
vibrations of rotors on underdamped flexible foundations. By including the energy loss 
associated with damping factors checked by experiment, they show close agreement be- 
tween theory and observation. Gunter’s work on a flexible-foundation analysis in refer- 
ence 3 reveals that proper foundation design can provide considerable improvement in 
rotor stability characteristics. 

The present report extends the methods of reference 1 to include vibratory motion of 
the foundation. It derives the general frequency equation, presents the results of a para- 
metric analysis of rotor dynamic frequencies and critical speeds, and gives some se- 
lected examples. The purpose of this report is to aid in preliminary design of turboma- 
chinery by providing guidance to the location of potential synchronous and nonsynchronous 
critical speeds. 

The work presented here covers a wide range of rotor shapes, rotor speeds, and 
foundation shapes. The centers of gravity of the rotor and foundation in the present 
model remain fixed midway between the bearings. Damping and unbalance forces are 
neglected in this analysis to isolate the effects of shape, speed, and spring- constant var- 
iations. 


SYMBOLS 


A, B, C, D, E, F, G, H 

a l> a 2> a 3’ a 4 

I 

■i 

kj/2 

k 2 /2 

l 

M 

M* 

P 

r 


S 

t 


vibration amplitude, m 

constant functions of 7 Tq, tt§, and i T q 

2 

polar moment-of-inertia of rotor, kg-m 

2 

diametral moment-of-inertia of rotor, kg-m 

2 

moment of inertia of foundation, kg-m 

o 

linear spring constant of rotor bearings, kg/sec 

2 

linear spring constant of foundation, kg/sec 

distance between bearings, and also between foundation mounts, m 
rotor mass, kg 
foundation mass, kg 
rotor precession frequency, rad/sec 
critical- speed ratio, r = = p/w 

square of rotor frequency parameter, 
time, sec 



2 



Yl’ Y 2 lateral distance from arbitrary reference frame to rotor bearings, m 

y'p y'g lateral distance from arbitrary reference frame to foundation mounts, m 

ZjZg lateral distance from arbitrary reference frame to rotor bearings (perpendicu- 

lar to y^ and y 2 ), m 

z’jzlj lateral distance from arbitrary reference frame to foundation mounts (perpen- 
dicular to y’j and y^), m 

77 ^ frequency parameter, p^/|/k^/M = Vs 

7^2 rotor disk effect, I^/MZ 2 

TTg rotor polar-to-diametral moment- of '-inertia ratio, I/I^ 

7 rotor rotational- speed parameter, u>y|/k^/M 
7Tg foundation-to-rotor moment-of-inertia ratio, Ij/Ij 

77g foundation-to-rotor spring- constant ratio, kgA^ 

77 ^ foundation-to-rotor mass ratio, M’/M 

a) rotor rotational speed, rad/sec 

Subscript: 

cr major critical 

Superscripts: 

first derivative with respect to time 
second derivative with respect to time 


DESCRIPTION OF ROTOR -FOUNDATION MODEL 

Figure 1 is a schematic representation of a rigid rotor with flexible bearings on 
pedestals attached to a flexible foundation. This study considers motions in the lateral 
mutually perpendicular directions y, y' and z, z' directions. Consideration of axial 
motion is excluded. 

No damping or unbalance forces are assumed in the model of figure 1. The centers 
of gravity of the rotor and foundation masses are assumed to be midway between the ro- 
tor bearings. All bearings are assumed to be linear springs. In both lateral directions 
the rotor- bearing spring constant is kj/2 and that of the foundation is kg/2 (see fig. 1) 

To this extent, there is one less variable in the present analysis than in that of ref- 
erence 1, because the center-of-gravity location is fixed. An additional variable is 
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Figure 1. - Schematic sketch of a rigid rotor on a flexible foundation. 

added, however, due to foundation flexibility in terms of the spring constant kg. Fur- 
thermore, the presence of the vibrating foundation mass adds two more variables - its 
mass M' and moment of inertia 1^. The net result is that the present analysis is com- 
plicated by two more variables than that of reference 1. 

Rotor slenderness ratio is introduced into the mathematical analysis in terms of the 
rotor polar moment of inertia I and the diametral moment of inertia Ij . The rotor 
moments of inertia are expressed as the following nondimensional parameters: 

Rotor disk effect: 



Rotor polar-to-diametral moment-of-inertia ratio: 



( 1 ) 


( 2 ) 


The parameters and ir^ are small for pencil shapes, and larger for disks. The 
parameter jjg is the disk effect used in reference 4. For a concentrated mass, = 0, 
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and for a disk having all its mass distributed over a large radius, 7 r 2 = °°. The ratio 
7 Tg = 0. 0149 for a rod whose length is 10 times its diameter. For a disk of negligible 
thickness, tt ^ = 2 . Therefore, configurations become more like disks as 7 ^ and 713 
increase. 

ANALYSIS 

This section presents the equations of motion for the simplified model shown in fig- 
ure 1. Solution of the equations of motion leads to a general frequency equation. This 
equation is applied to various special conditions of geometry and speed to cover the wide 
range of conditions used in the examples. 

The object of this analysis is to study the frequency and critical- speed characteris- 
tics of the model in figure 1. Therefore, no attempt is made to solve the equations for 
vibration amplitudes. 

Equations of Motion 

Derivation of the following equations of motion was motivated by the work of 
Timoshenko in reference 5. The equations include the assumption of small rotor and 
foundation displacements: 

Rotor: 

Myj + My 2 + k 1 (y 1 - y^) + kj(y 2 - y ? 2 ) =0 ( 3 ) 

Mzj + Mz‘ 2 + kjCzj - Zj) + kj(z 2 - Zg) = 0 (4) 

k k 

Icu(y 2 “ yj) ~ -^(^2 ” z *l^ “ — ( z 2 ~ z 2^ ^ — ( z j ™ z l) ~ ® (®) 


k ^ k 

Ico(z 2 - Zj) + I x (y 2 - y x ) + — — (y 2 - y' 2 ) - — — (y A * y^) = 0 (6) 

,w 4 4 

Foundation: 

(7) 


M’y’j + M’y ’ 2 - k 1 (y J + y ? ) + (kj + kgXy’j + y^) = 0 
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( 8 ) 


M'z'i + M’z’g - k 1 (z 1 + z 2 ) + (kj + k 2 )(z , 1 + Zg) = 0 


(k^ + k 0 )r 

-tyt'2 - * 1 ) + - 1 2 


k- l l 


4 t*i - z y - -7- - z 2> = ° 


(9) 


, (k, + k «,)£ 2 k ^ 2 

r i( y' 2 - fl ) - 1 2 (y\ - y’j) ♦ 4- (y, - y 2 ) - o 


( 10 ) 


Frequency Solutions 


Assumed solutions . - Den Hartog (ref. 4) solves problems of two masses joined by 
springs by assuming that both masses execute harmonic motions with the same fre- 
quency. This assumption is conservative, because it represents a resonant condition in 
which the motions of the two masses reinforce each other. This is the most critical 
condition, and is the one to be avoided if possible. 

Accordingly, the same frequency p appears in all of the following assumed solu- 
tions: 


n 




z i 

z 2 


A 


*2 

z i 


= A sin pt 

(11) 

= B sin pt 

(12) 

= C cos pt 

(13) 

= D cos pt 

(14) 

= E sin pt 

(15) 

= F sin pt 

(16) 

= G cos pt 

(17) 

= H cos pt 

(18) 
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Coefficient matrix . - When equations (11) to (18) are inserted into equations (3) 
to (10), the following set results: 


A(kj - Mp 2 ) + B(kj - Mp 2 ) - Ekj - Fkj = 0 


( 19 ) 


C(kj - Mp 2 ) + D(kj - Mp 2 ) - Gk 1 - Hkj = 0 


( 20 ) 


-A(Iwp) + B(Iwp) + c(^- - Ijp 2 ^ - " !iP 2 



= 0 (21) 


-a(^_ - ^p 2 ^ + " x iP 2 ) + C(lwp) - D(Iwp) + E 



= 0 (22) 


-Ak 1 - Bk x + E [(kj + k 2 ) - M'p 2 ] + F[(k x + k 2 ) - M’p 2 ] = 0 (23) 


-Ck 1 - Dkj + G[(kj + k 2 ) - M’p 2 ] + H [(kj + kg) - M’p 2 ] = 0 (24) 



+ G 


(kj + k 2 )r 2 

-T 2 - - V 


- H 


(k 1+ k 2 )i 2 

ljp 


= 0 (25) 



- E 


+ k 0 )l* 9 

-i_ . lip 2 


+ F 


(k< + kn)l 2 

-M i\p 2 


= 0 (26) 


In matrix form, this set is 
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Determinant reduction . - Nontrivial solutions of A, B, C, D, E, F, G, and H in 
equations (19) to (26) require that the determinant of the coefficient matrix given in equa- 
tion (27) vanish. Reduction of this determinant leads to the general frequency equation 




Equation (28) is of the 16^ order in p. There must therefore be 16 solutions, al- 
though not all are unique. In equation (28), the determinant and its multiplying factor 
each contains eight solutions. 

Nondimensional parameters. - Expressing the variables in equation (28) in nondi- 
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mensional form is especially helpful in solving the equation and presenting the results. 
The nondi mensional parameters and their nomenclatures are defined as follows: 

Precession frequency parameter: 


n l = 



Rotor disk effect: 


7T n 


Mr 


Rotor polar-to-diametral moment-of-inertia ratio: 


*3 = 


Rotor-rotational- speed parameter: 


H = 



Foundation-to- rotor moment-of-inertia ratio: 



Foundation-to- rotor spring- constant ratio: 



(29) 


(1) 


( 2 ) 


(30) 


(31) 


(32) 


Foundation-to-rotor mass ratio: 
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(33) 



An objective of this analysis henceforth will be to identify frequency solutions for 
7 Tj instead of p. It will be seen that definition of another parameter S simplifies the 
appearance of high-degree polynomials by halving their orders: 



Nondimensional general frequency equation . - With the use of the nondimensional 
parameters defined by equations (1), (2), and (29) to (33), the general frequency equation 
(eq. (28)) assumes the following nondimensional form: 



Mode Frequency Equations 

In equation (35), both the determinant and its multiplying factor contain valid solu- 
tions. Solutions for obtained from this factor represent the bouncing-mode natural 
frequency. The rotor rotational-speed parameter it 4 and moment-of-inertia ratios 7 ^, 
TTg, and 7 Tg do not appear in this factor. Therefore, the bouncing mode solutions are 
the natural frequencies that exist whether the rotor is rotating or not. 

The determinant in equation (35) yields the conical- mode solutions, which are func- 
tions of rotor speed and the moment-of-inertia ratios. The foundation-to- rotor mass 
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ratio 7 Trj, however, is absent from the determinant. Masses M and M' exist in the 
determinant, of course, in terms of the moment- of-inertia ratios. The conical mode 
consists of the rotor precessions performed as the rotor axis describes the surface of 
an imaginary cone. The study of the conical- mode solutions is a primary concern of this 
analysis. 

Bouncing mode . - Equating to zero the factor outside the determinant in equation (35) 
and using equation (34) yield the bouncing-mode frequency equation 

^ S 2.(l + iJ S+ l = 0 (36) 

w V ^6 v 

The solution of this quadratic expression is 



Values of S greater than one result from using the positive sign in this expression. The 
negative sign yields values of S less than or equal to one. 

Equation (37), by way of equation (34), represents 4 of the total of 16 solutions for 
p or contained in equations (28) or (35). Because the factor outside the determinant 
in equation (35) is of second degree, there exists a duplicate set of the four solutions for 
7Tj contained in equation (37). Of these eight solutions, only two are unique. 

Conical mode . - The conical- mode frequency equation is obtained by equating the 
determinant in equation (35) to zero. The result is 
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Solving this quartic equation for S is cumbersome. It is more convenient to solve 
equation (38) for the shape-speed parameter ir^. The resulting expression is 




+ 



S 
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Further simplification is possible by multiplying equation (39) by and considering 

(Sffg) and (^ 4 ^* 2 ) as the variables. Therefore, equation (39) becomes 



The arrangement (Sn^) in effect reduces the number of variables by 1. 
As (Sjrg) -0 in equation (40), the relation 
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( 41 ) 



exists with the aid of equation (34). This relation is approached asymptotically by back- 
ward precession in the low-frequency set. 


As (S7r 2 ) 



in equation (40), the ratio 



= 1 


(42) 


is approached for all shapes, irrespective of and 7Tq. The relation expressed by 
equation (42) applies to the condition of forward precession in the high-frequency set. 
When the spring- constant ratio jtq — 0, an approximate form of equation (40) is 



(Sir 2 ) = (Sir 2 ) 


2 




(43) 


To provide physical significance of the zero condition for ttq, it may be considered that 
it is approached by rotors in aircraft or spacecraft (floating foundation). Very hard 
spring rotor -bearing mounts also approach this condition. 

A remarkable simplification to the solution of equation (43) is available by noting 
that the right side is a perfect square. A convenient form is 


where 

and 



5 = 1 for forward precession ^ 

► 

6 = -1 for backward precession. 


(44) 


(45) 
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Figure 2. - Coefficient for abbreviated solution (for use in eq. (44)). 


The technique of displaying the frequency equation in the form of equation (44) was sug- 
gested by Den Hartog in reference 4. It was also used in reference 1. 

Figure 2 presents the variati on with the foundation-to-rotor moment-of-inertia ratio 
77 ^ of the coefficient ^1 + (l/ffg) in equation (44). The sensitivity becomes extreme 
when 7 rg < 1, which represents a large rotor on a small foundation. 

When the spring- constant ratio tt 6 and square of the adjusted frequency parameter 
S both become large (high-frequency set), equation (40) reduces to the relation 


or 


S7t 2 = — (46) 

2 H 



( 47 ) 
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which is a parabola in 77^ V"2 and 77g. A large value of ttq approaches the condition 
for a firm foundation. This condition also applies to very soft rotor-bearing mounts on 
a flexible or firm foundation. The expression in equation (47), which will be discussed 
on page 29 in connection with figure 7, applies to the high-frequency set of solutions. 
Equations (46) and (47) also result from equation (40) when the adjusted shape-speed 

parameter 7 13 
When TTg 
tion (40) becomes 


1 — \ 


( 77-4 y 77^ J becomes large, as by high rotor rotational speed, 
aecomes large with S77- 2 remaining finite (low-frequency set), equa- 


7T 


2 

3 



(Stt 2 ) = (Stt 2 ) 2 



(48) 


It is interesting to note that equation (43) reduces to the expression in equation (48) as 
7Tg — °°. This condition pertains to small rotors on large foundations. 

Again, because of the perfect square on the right side of equation (48), its solution 
is 


677g2 




(49) 


This result is identical to that reported in reference 1 for the case with the center of 
gravity midway between the bearings. 

At zero rotor rotational speed, equation (44) reduces to 


*1 



(50) 


and equation (49) becomes 





1 

2 


(51) 


These solutions are entirely apart from the bouncing-mode frequencies of equation (37) 
because of the fundamental difference in their derivations. Certain geometrical condi- 
tions do exist, of course, in which the frequencies of the two modes could coincide. The 
result of equating expressions for S obtained from equations (37) and (50) is 
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Forcing the zero-speed and bouncing-mode frequencies to be equal in this way is pri- 
marily an artificial procedure. It is hoped that pointing out the artificiality will help 
avoid confusion between the zero-speed and bouncing-mode frequencies. 


Critical Speeds 

A critical speed is defined as any rotor rotational speed at which vibration amplifi- 
cation occurs. The bouncing mode is excluded from this classification because it is in- 
dependent of rotor speed. The critical- speed ratio is 


r = 


^1 = p 

7J4 w 


(53) 


When this relation together with equation (34) is used in equation (40), the following 
fourth- order equation in Sit 2 results: 


256 


'V?Al 

^6/ \r 2 


- l] (Sjt 2 ) - 128 


1 + 


_5 _ JL 
f 6 *6 


(S7r 0 r 


+ 16 


1 + i 


- 1 - 


^5 _ f*5 _ ^5 _ _2_ _ ± 

rr R 2 2 7r„ 2 

67 r 6 *6 6 ^6 


(Sn o y 


+ 8{l + 2+X|(S7r 9 ) -1 = 0 

n 6 % 


(54) 


A concise form of this expression is 


a 4 (S7r 2 ) 4 + a 3 (S7r 2 ) 3 + a 2 (S7r 2 ) 2 + a^S^) -1 = 0 


(55) 
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When a.p a.^, a^, and a^ are constant, equation (55) yields four solutions of the type 

Sn 2 = K 2 (56) 

or 

(57) 

The eight roots given by equation (57) are all real when ir^/r < 1, as equation (54) 
shows. Six roots are real and there is a complex conjugate pair when v^/r 1. The 
six real roots comprise two sets of mirror images. 

The mathematical condition for major critical speed is obtained from equation (53) 
by setting 


r = — = 1 (58) 

*4 

When this relation is used in equation (54), the result is 



As 7Tg — 0, the solution of equation (43) under the condition of equation (58) is 
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At the other extremity, when 77g = equations (48) and (58) yield 

(it a cr f~ 2 ) = - - l / — ( 61 ! 

\ 4, cr ? 2 /^gsoo 2 ^ 1 - 6 it 3 

Equation (60) reduces to equation (61) when 77^ becomes infinite. 

RESULTS AND DISCUSSION 

As previously mentioned, two classes of solution to the general frequency equations 
(eqs. (28) and (35)) exist. The determinant is a function of rotor rotational speed and 
shape; whereas its multiplying factor is independent of these parameters. This section 
presents graphical solutions and discussions of both the bouncing and conical modes of 
the frequency. 


Bouncing Mode 

The bouncing - m ode f r equency , obtained from equation (37), is the counterpart of the 
trivial solution p/ y^k^/M = 1, which is one of the solutions of the frequency equation 
derived in the firm-foundation analysis of reference 1 when the center of gravity is mid- 
way between the bearings. 

Figure 3 is a logarithmic plot of equation (37) that covers a range of the spring- 
constant ratio 770 essentially from zero to infinity. The ordinate is the dimensionless 
frequency parameter 77-^ = /M defined by equation (29). The mass ratio 

TTrj = M’/M varies over a range from one-half to infinity in figure 3. The upper set of 
solutions shown in figure 3 results from using the plus sign in equation (37). The minus 
sign yields the lower set. 

The absence of rotor rotational speed and moment of inertia in figure 3 is worth 
discussing. Each set of two points on the curves of this figure at a given value of ttq 
represents the bouncing- mode natural frequencies of one particular rotor configuration. 
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Figure 3. - Bouncing-mode natural-frequency plot (eq. (37)). 



The magnitudes and existence of these frequencies are completely dissociated from the 
stationary or rotating condition of this particular rotor. Only the spring- constant ratio 
7Tg and mass ratio 7 affect these solutions, as equation (37) and figure 3 both show. 

As the spring- constant ratio -itq becomes infinite, figure 3 and equation (37) with 
the minus sign reveal that the bouncing- mode frequency parameter approaches one. 
Likewise, as the mass ratio 7 — ■«, figure 3 and equation (37) with the plus sign show 
that ttj becomes one. Physically, these two conditions correspond to the firm- 
foundation result 77 ^ = 1 found in reference 1 . 

At the other extremity in figure 3, where ttq — 0, the lower set of solutions van- 
ishes. The upper set approaches the value 1 + (l/jr»), as equation (37) shows. 


Conical Mode 


The main discussion in this report is concerned with the conical- mode natural fre- 
quency, which is obtained mathematically when the determinant in equation (35) vanishes. 



=3 

w 


g 

1 0 
=5 
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■§ 
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xa 
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Figure 4. - Abbreviated solution (eq. (44) for 7r 6 = 0; 
eq. (49) for' 7T 6 = oo). 
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The solution is given by equation (40). It is evident that the foundation-to-rotor mass 
ratio 7 Trj is absent from this solution. 

Prior to a discussion of the graphical solution of equation (40), consideration of 
special cases of the spring- constant ratio {ttq = 0 and ttq = °°) is instructive. Applica- 
tion to aircraft and spacecraft for the condition ttq = 0 was cited in connection with the 
derivation of equation (44) (see p. 13). The case tIq - 00 represents the firm- foundation 
model of reference 1 . 

Abbreviated solution . - Equation (44) is plotted in figure 4 with the factor 


2/Vl + (1A 5 ) multiplying ^ in the ordinate and multiplying 77 3 (rr^ 
abscissa. Inspection of equation (49) shows that figure 4 is also a graphical solution for 
this equation if the foundation-to-rotor moment-of-inertia ratio 7 Tg is assumed infinite 
in the coordinates. Therefore, a feature of figure 4 is that it applies to each of the two 
extreme cases (floating foundation, 77 g = 0 and firm foundation, ttq- °°). The set of 
curves in figure 4 is identical to a set in reference 1, with appropriate adjustment of 
coordinates. 

The advantage of the form of solution given by figure 4 is that all geometries can be 
displayed by a single two-branch curve. The upper branch shown by the solid curve 
represents forward precession. The dotted curve represents backward precession. 
Figure 4 Shows that as the adjusted shape- speed parameter increases, as by rising rotor 
rotational speed, the forward-precession frequency increases and the backward fre- 
quency decreases. 

The point of intersection of the two branches in figure 4 locates graphically the 
zero-speed condition given analytically by equations (50) and (51). These results have 
therefore shown that a stationary rotor on a floating or firm foundation has two natural 
frequencies - one bouncing- mode frequency from equation (37) and one conical- mode fre- 
quency from equations (50) or (51). 

From a comparison of equations (37 ) and (52), it is apparent that figure 3 also pre- 
sents the geometric constraints for ^ ^1 + j ^2 that would be required for the 


zero-speed and bouncing-mode frequencies to be equal. 

General solution . - Some example solutions of the general frequency equation (40) 
are presented graphically in figure 5. Representative samples appear in figure 5 with 


adjusted frequency parameter y tj^ and adjusted shape- speed parameter 773 ^774 
as coordinates. Presence of the disk effect 77g in the coordinates reduces by two- 
the number of plots needed to show the amount of information in figure 5 The 
foundation-to- rotor moment-of-inertia ratio 775 = Ij/I^ is assigned values of 0. 5, 1, 
and 3 in each part of figure 5. The spring- constant ratio 7 Tg = kgA^ assumes constant 
values of 1/3, 1, and 2 in figures 5(a), (b), and (c), respectively. 

For every value of ti^ on each plot in figure 5, there is a high-frequency set of so- 


til 

■thirds 




(a) Foundation -to -rotor spring -constant ratio, 
tt 6 = (kj/kj) - 0.3333. 

Figure 5. - General solution plots (eq. (40)). 


lutions and a low-frequency set. For firm foundations, reference 1 reports only one set 
when the center of gravity is midway between the bearings. For given geometry, fig- 

except when 

curves represent forward precession, while the dotted curves represent backward pre- 
cession. The backward-precession curves rightly belong in the fourth quadrant but ap- 
pear in the first quadrant for convenience. A mirror image of a general solution plot 
also exists in the second quadrant, but only repeats the information given in figure 5. 
These eight roots from figure 5 and its mirror image comprise the total implied by 


ure 5 shows four distinct solutions for 


'iV** 




= 0. The solid 
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(b) Foundation -to-rotor spring -constant ratio, 
w 6 = (k 2 /k x ) = 1. 

Figured, -Continued. 


equation (38). Because of the presence of 7 ^ and 7 in the coordinates, each set of 
four points in figure 5 at a particular rotor rotational speed represents a single rotor 
configuration. 

Comparison of figures 4 and 5 reveals a primary result of this investigation. This 
is an assessment of the extent of the error made by approximating a flexible-foundation 
configuration by a firm-foundation model. It is evident immediately that the flexible- 
foundation model predicts the existence of four frequency solutions, whereas the firm- 
foundation model predicts only two. 
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Adjusted shape-speed parameter, ^3^4 $ 2 )" 



(c) Foundation-to-rotor spring-constant ratio, 

tt 6 = (kg/kj) = 2. 

Figure 5 . - Concluded. 


Table I compares quantitative results predicted by the two models. This table shows 
the approximate agreement between the flexible-foundation results and the two corre- 
sponding solutions of the firm-foundation model. Furthermore, such agreement exists 
over a wide range of spring- constant ratio 7Tg, as will be shown by the discussion on 
page 29. 

Figure 5 illustrates the following trends when the effects of ir^, Vgj 71 an d Vg 

are each isolated from all the others. 

(1) When the rotational speed of a given rotor rises, the forward precession fre- 
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TABLE 1. - COMPARISON OF FREQUENCY PARAMETERS PREDICTED 
BY FLEXIBLE- AND FIRM- FOUNDATION MODELS 


Foundation 

model 

Fig- 

ure 

*5 

H 


High- 

frequency set 

Low- 

frequency set 







\ f ~2 


\^2 

Flexible 

5(c) 

1 

2 

2. 5 

2.6 

0.88 

0.84 

-0.065 

Firm 

4 

OO 

OO 

2.5 

2.6 


— 

-.10 


a Forward precession. 
^Backward precession. 


quency increases and the backward frequency decreases. 

(2) The frequency parameter 7 Tj is inversely proportional to the square root of the 
disk effect y 

(3) The forward-precession frequency increases and the backward frequency de- 
creases as the rotor polar-to-diametral moment-of-inertia ratio rises. 

(4) The precession frequency decreases as the foundation-to-rotor moment-of- 
inertia ratio Fg increases. 

(5) When the foundation-to-rotor spring- constant ratio ttq increases, the preces- 
sion frequency increases. 

When the adjusted shape-speed parameter becomes large, all the 

forward-precession high-frequency curves in figure 5 approach as an asymptote the re- 
lation expressed by equation (42). Figure 5 also shows the asymptote approached by 
backward precession in the low-frequency set, as expressed in equation (41). 

In figure 5, there is a zero-speed condition for each of the two sets, given by the 
intersections of the branches. Therefore, a stationary rotor on a flexible foundation 1ms 
four natural frequencies, including the two bouncing- mode solutions. 

Frequency plots. - A graph of the variation in adjusted frequency parameter V*2 
with adjusted rotational- speed parameter ir^ is called a frequency plot. It is easily 
constructed from a general solution plot by specifying values for the rotor moment-of- 
inertia ratio, jig for use as a parameter as in figure 6, In this figure, which was de- 
rived from figure 5(c), Vg is 3 and tTq is 2. Therefore, figure 6 is less general and 
presents less information than figure 5(c). 
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Precession 



Figure 7. - Spring-constant plots. 



Shape-speed parameter 



(b) 7T 2 = I 2 /M^= I; 7T 5 « Ii/I 3 
Figure 7. - Concluded. 







Frequency plots are convenient for locating critical speeds by intersections of 
critical- speed ratio curves and frequency curves. Several critical- speed ratio curves 
appear in figure 6. Their significance will be treated in the nonsynchronous critical- 
speed discussions (see pp. 39 and 40). Each constant 77 g curve on a frequency plot 
represents the behavior of a particular rotor, and shows its passage through various 
critical speeds during rotor rotational- speed changes. 

The gradual slope of the low- value 77g curves in figure 6 shows that variation in 
rotor rotational speed has little effect on the precession of pencil- shape rotors. This ob- 
servation is also true for high- and low-frequency backward precession of all shapes. 

Spring- constant plots . - An overall perspective of the effect of variations in the 
foundation- to- rotor spring- constant ratio Vg on the frequency parameter is provided 
by spring- constant plots, two of which appear in the logarithmic graphs of figure 7. This 
figure uses the shape- speed parameter as the parametric variable. Figure 7(a) 

applies to any geometry in which 77 g = 0. 05 and 77 g - 0. 5, while 7(b) applies to any 
geometry in which n g = 1 and 77 g = 1. Figure 7 covers a range in spring- constant ratio 
77 g from essentially zero to infinity. Figure 4 provides the totality of solutions for both 
the lower and upper boundaries of 77g in figure 7. 

The flatness of the curves in figure 7 account for the agreement between the firm- 
and flexible-foundation models shown in table I. Figure 7 shows close agreement of the 
high-frequency forward-precession curves over a wide range of spring- constant ratio 77g 
with the low-frequency forward- precession curves at their firm- foundation limits. Like- 
wise, the low-frequency backward-precession curves in figure 7 are insensitive to varia- 
tion in 77 g above values for 77g of about 1. 

As the spring- constant ratio 77 g increases, the curves of constant 77 ^ 77 ^ of the high- 
frequency set approach asymptotically the relation given by equation (47), as shown in 
figure 7. On a logarithmic plot of 77 ^ against 77g, equation (47) maps as a straight line. 

As 77 g increases, the curves of the low-frequency set in figure 7 approach values 
given by figure 4 and equation (49) with zero slope. 

Thus, for floating foundations, the solutions provided by figure 4 correspond to the 
high-frequency set. For firm foundations, the figure 4 solutions correspond to the low- 
frequency set. 

Figure 8 is an expanded portion of figure 7(a), showing the region about 77 q = 0 in 
rectangular coordinates. 

To summarize, plots such as those in figure 7 provide solutions for the entire range 
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Shape-speed parameter, 



Figure 8. - Enlarged section of spring-constant plot of fig- 
ure 7(a). 7 r 2 = I iim? - 0. 05; tt 5 = iyij = 0. 5. 


of geometry, speed, and spring- constant ratio. Use of only three values of spring- 
constant ratio 77 g in calculation of the general frequency equation (40) suffice in produc- 
tion of any one of these plots for the following reasons. The Den Hartog-type plot of fig- 
ure 4 establishes the lower and upper extremities of itq in figure 7. The high-frequency 
curves approach asymptotically the expression given by equation (47) as ttq increases. 

Asymptotic limit . - Figure 9 is a graphical representation of the asymptotic limit 
given by equation (47). This figure shows that the parabolic relation between 7 and 
?Tg is linear with a slope of 0. 5 for all values of Trgflg on a logarithmic plot, as taking 
the log of equation (47) verifies. As is evident from figure 7, the asymptotic limit forms 
the backbone of the spring- constant plots. 
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Combined 
moment-of-inertia 
parameter, 
7^775 = I j /A/ll 2 



Figure 9. - Asymptotic limits (eq. (47)). 

Major Critical Speeds 

A practical application of the theory and maps developed in this report is to use 
them to locate major critical speeds. Knowledge of major critical- speed conditions is 
important because large vibration- amplitudes occur at major critical speeds in the ab- 
sence of damping. Major critical speed is defined by equation (58) 

Major critical- speed conditions are most easily visualized by referring to a fre- 
quency plot, such as figure 6. On a frequency plot, intersections of the critical- speed 
ratio curve r = 1 with the frequency curves establish the major critical speeds. The 
trace of a major critical-speed curve on a general solution plot, such as figure 5, is 
given by 



Effect of shape . - A less pictorial but more mathematically rigorous approach to 
locating the major critical speeds is provided by equation (59). Graphical solutions of 

this equation appear in figure 10, which uses ff 4, cr V^2 anc * ^3 as coor dinates. The 
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Rotor polar-to-diametral moment Df-inertia ratio, ** (I/Ij) 

(a) Foundation-to-rotor spring-constant ratio, 

7r 6 = (k 2 /k i) - 0.3333. 

Figure 10. - Adjusted major critical -speed plots (eq. (59)). 


parameter varies in figure 10 to show the effect of rotor shape on major crit- 
ical speed. The spring- constant ratio is fixed in figure 10(a) at 0. 333, in 10(b) at 1.0, 
and in 10(c) at 2. 0. 

Trends (2) to (5) enumerated in the General solution section (pp. 24 and 25) also ap- 
ply to major critical speeds, as figure 10 shows. 

All the forward-precession curves of the high-frequency sets in figure 10 approach 
infinity as ir 3 — 1. This circumstance is related to the discussion of the roots of equa- 
tion (54) (see p. 17). Therefore, because no real solutions exist for S: 1, forward 
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(b) Foundation -to-rotor spring -constant ratio, 
7r 6 = (k 2 /k 1 ) = 1. 


Figure 10. -Continued. 


major critical speed of the high-frequency set can be avoided by using such designs. If 
this is not feasible, designing with 7 as close to one as possible may yield a major 
critical speed well above the design rotor rotational speed. 

Selective geometrical design can not so easily avoid backward major critical-speed 
resonances, as shown by the fact that the dashed curves in figure 10 remain finite over 
the entire range of ijg. The presence of ^71^ in the ordinate of figure 10, however, 
suggests that if a design with a low value of 7^ is used (pencil shapes), both the forward 
and backward major critical speeds may be high enough to exceed design speed. 

Effect of spring- constant ratio . - Figure 11 presents the effect of spring- constant 

ratio ?7g on ir^ cr ^^2 w ^h curves of constant values of 77 5. A pencil- shape rotor 
77g = 0. 1 appears in figure 11(a), and a disk fig = 2 is represented in figure 11(b). 

Figure 11 gives an overall perspective of major critical- speed conditions. The 
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(cl Foundation-to-rofor spring-constant ratio, « (k 2 /k 1 ) = 2. 


Figure 10. - Concluded. 


curves at the lower extremities of ttq are approximated by equation (60), and at the 
upper by equation (61). The asymptotes approached by the high-frequency set are plotted 
from equation (47). 

No high-frequency forward-precession curves appear in figure 11(b) because r < ir^ 
(seep. 17). As tTq becomes large, the low-frequency curves in figure 11(a) converge 
into one forward- and one backward-precession curve (eq. (61)). The backward- 
precession curves of figure 11(b) behave like those in figure 11(a). The low-frequency 
forward-precession curves of figure 11(b) merge into their respective asymptotes as 7 Tg 
becomes large, as equation (47) specifies. 

There are no regions in the quadrant represented by figure 11 in which critical- 
speed solutions are forbidden. In figure 11(a), as — 0 the low-frequency curves ap- 
proach the constant value of equation (61) for all values of ttq. This observation also is 
true for the low-frequency backward- precession curves of figure 11(b). But when 
TTg > r (fig. 11(b)), the low-frequency forward-precession curves approach the ordinate 
axis (not shown) as n§ — 0. 
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Figure 11. - Effect of spring-constant ratio on major critical speed. 








As Tig increases, the low-frequency set covers the lower region of figure 11 all the 
way down to the abscissa axis (not shown). As Vg — 0, curves of the high-frequency set 
approach infinity in figure 11 . At the other extreme, they approach the constant value 
given by equation (61) as Vg becomes infinite. 

Because of the trends observed in figure 7, it was expected that for a given geom- 
etry the major critical speeds shown in figure 11 would increase with an increase in irg. 

Rising rotor rotational speed is represented in figure 11 by a traverse upward on 
the map at constant ttq. Such a movement from a Vg curve of the low-frequency set to 
the corresponding Vg curve of the high-frequency set locates successive major critical 
speeds. 

Abbreviated major critical- speed plot . - Just as the abbreviated frequency plot in 
figure 4 set the extreme boundaries of figure 7, the abbreviated major critical - speed 
plot in figure 12 does likewise for figure 11. With the ordinate 2ir^ .crV^/A + (1A 5 ), 
figure 12 is a plot of equation (60) for which Uq = 0 (floating foundation). When 7Tg is 
set at infinity, figure 12 represents equation (61), for which itq = °° (firm foundation). 

The rotor moment- of-inertia ratio is the abscissa coordinate in figure 12. It varies 


^6 = G n 6 = 00 



Figure 12. - Abbreviated major critical -speed 
plot (eq. (60) for 7r 6 * 0; eq. (61) for 7^ = °°). 
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TABLE n. - COMPARISON OF MAJOR CRITICAL SPEEDS PREDICTED BY 


FLEXIBLE- AND FIRM -FOUNDATION MODELS 


Foundation 

model 

Shape 

Figure 

*5 

*6 

ff 3 

High- 

frequency set 

Low- 

frequency set 

n 4, cr V^2 

v 4, cr V^2 

%, cr V^2 

\ cr V*2 

Flexible 

Pencil 

10(c) 

1 

2 

0.1 

0.93 

-0.91 

0.40 

-0.37 

Firm 

Pencil 

12 

oO 

oo 

.1 

(1) 

(-D 

.52 

-.47 

Flexible 

Disk 

10(c) 

1 

2 

2 

— - 

-.88 

,88 

-.23 

Firm 

Disk 

12 

OO 

OO 

2 


(-1) . 

— - 

-.29 


a Forward precession. 
°Backward precession. 


from zero to the value for a disk (2). Again, as in figure 10, the high-frequency curve 
approaches infinity as — 1, beyond which this curve does not exist. 

Comparison with firm-foundation results . - Table II presents a comparison of the 
major critical speeds predicted by the two models for both pencil-shape rotors and disks. 
In table n, the two solutions not predicted by the firm-foundation model are approximated 
by using plus and minus the firm-foundation bouncing-mode solution. The reasonably 
good agreement between the two sets of results for both pencil-shape rotors and disks 
implies good agreement for all intermediate shapes. Table II has thus shown that the 
firm-foundation results from figure 12 with the addition of this model's bouncing-mode 
solutions (±1) provide an easily obtained but good approximation to flexible-foundation re- 
sults. 


Nonsynchronous Critical Speeds 

In references 6 to 9, Yamamoto reports on experimental work showing that large 
amplitude increases occur at rotor rotational speeds other than the major critical. His 
work with rotors supported by ball bearings revealed that bearing defects cause nonsyn- 
chronous precession. He found two classes of nonsynchronous precession associated 
with double-row ball bearings (refs. 8 and 7), and two other classes associated with 
single-row types (refs. 8 and 9). 

With double-row ball bearings, nonsynchronous critical speeds occur both above and 
below the major critical speed. All nonsynchronous critical speeds associated with 
single-row ball bearings occur at rotor rotational speeds greater than the major critical. 
The spring characteristics of double- row ball bearings are linear, and those of the 
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single-row type are nonlinear. 

Double-row ball bearings . - Nonuniformity in ball diameters is a cause of nonsyn- 
chronous motion at rotor rotational speeds greater than the major critical. Yamamoto's 
work disclosed that the most serious conditions of this type occur when 

r = — — (63) 

2. 65 

for forward precession, and 

r = - — (64) 

4.1 


for backward precession. 

Frequency plots, such as the example shown in figure 6, are ideally suited for dis- 
playing critical-speed conditions. Curves of constant critical- speed ratio r could also 
be shown on the general solution plots of figure 5. 

Curves representing equations (63) and (64) appear in figure 6. Only the intersec- 
tions of equation (63) with the forward-precession curves are meaningful. Similarly, the 
intersections of equation (64) with backward-precession curves are meaningful. Dots in 
figure 6 denote valid intersections. For a given rotor, the nonsynchronous critical 
speeds specified by equations (63) and (64) are more than twice the major critical speeds. 

Yamamoto showed that noncircular inner and outer bearing races cause nonsynchro- 
nous motion at critical-speed ratios r of 2, 3, and 5. Because of the small amplifica- 
tion that occurs at r = 5, this condition is not serious. Curves depicting r = 2 and 
r = 3 appear in figure 6. According to Yamamoto, only the intersections with forward- 
precession curves are valid. Because all intersections occur well below the major crit- 
ical Speed, even low-speed rotors can experience this type of nonsynchronous precession 
during acceleration. 

Inspection of the curve for 7Tg = 2 in figure 6 reveals that no high-frequency 
forward-precession intersections occur for critical- speed ratios r < 2. This is a 
graphical illustration of the fact discussed in connection with equation (54) that one of the 
solutions for 7r^ Vi is imaginary when ir^/r > 1. Thus, as reference 1 reported, no 
high-frequency forward-precession critical speed can exist for geometries in which i 
equals or exceeds the critical- speed ratio. Low-frequency critical speeds exist, how- 
ever, for all rotor shapes. 

It is not likely that a set of double-row ball bearings will have only one of the two 
types of defects mentioned. It may therefore be presumed that nonsynchronous critical 
speeds defined by equations (63) and (64), and r = 2 and 3 will all occur. 
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Single-row ball bearings . - Yamamoto observed in reference 8 that nonsynchronous 
critical speeds of subharmonic oscillation appear only when the shafts are supported by 
single-row ball bearings. This motion, which appears only in forward precession, re- 
sults from nonsymmetry of the nonlinear spring characteristics of the bearings. 

Subharmonic motion is characterized by critical-speed ratios r = 1/2, 1/3, 

1/4, . . , 1/n, where n is a positive real integer. This phenomenon can occur only 
after the major critical speed is exceeded because r < 1 (see fig. 6). Also, because of 
the low value of r, pencil-shape rotors experience both high- and low-frequency motion 
of this type. Only low-frequency precession is possible in disks. Figure 6 shows curves 
of r = 1/2 and r = 1/3 to provide examples of this phenomenon. 

Yamamoto further observed in references 8 and 9 that when single-row ball bearings 
are used, two natural frequencies build up together at certain rotor rotational speeds. 

The absolute value of their sum or difference is related to rotor rotational speed by 



in which p. and Pj are any two of the four natural frequencies obtained from equa- 
tion (54), Yamamoto refers to this motion as summed and differential harmonic oscilla- 
tions. In reference 8, Yamamoto reports that vibration amplitudes from these sources 
may exceed those of the major critical speeds. From figure 5, it is evident that the 
condition i = j is trivial and applies only at zero rotational speed. 

Because there are four unique solutions p^, Pg, Pg, and p^ to the frequency equa- 
tion, there are 12 summed and differential combinations of the type specified in equa- 
tion (65). To display these combinations graphically, figure 13 presents the portion of 
the frequency plot of figure 6 pertaining to flg = 0. 4. Solutions p^ and Pg represent 
forward precession and Pg and p^ represent backward precession. The curve r = 1 
locates nonsynchronous critical speeds of this type in figure 13. Because the intersec- 
tions of all curves with the curve r = 1 in figure 13 are valid, it is unnecessary to use 
dots. 


At zero rotor rotational speed, solutions p^ and p^ coincide as do Pg and Pg. 
Therefore, there can be no intersections of the curve r = 1 with the curves | p ^ + p^ | 
and Ipg + Pgl, as figure 13 shows. These solutions are forbidden for any geometry 
from pencil shape to disk. Figure 13 shows that the maximum number of intersections 
is 10 for any shape. 

Yamamoto and Hayashi proved mathematically in reference 9 that only summed har- 
monic oscillations can occur if both frequencies pj and Pj represent forward preces- 
sion, or if both represent backward precession. Only differential harmonic oscillations 
can occur if one frequency represents forward and the other represents backward pre- 
cession. Table in reviews the combinations allowed and forbidden. 


40 



Adjusted fr 


Major 



Figure 13. - Frequency plot locating sum med-and- 
differential harmonic nonsynchronous critical 
speeds. 773 = I/Ij = 0.4; 7 r 5 = iyij * 3; 

H* k 2 /k l = 2 * 


TABLE IE. - ALLOWED AND FORBIDDEN SUMMED AND 
DIFFERENTIAL HARMONIC COMBINATIONS 


Allowed 

|pi + p 2 l 

bi -P3I 

bi - P4I 

|p 2 - P3I 

|p 2 - P4I 

|p 3 + p 4 l 

Forbidden 

Ipi -p 2 l 

|pi + p 3 l 

1 pi + P4 1 

|p 2 + P3I 

|p 2 + P4I 

IP3-P4I 




For disks, no |p^±Pj| intersections can exist, as figure 8 verifies. These inter- 
sections are nonexistent for any rotor shapes for which a 1. Consequently, for disks 
only 5 intersections of |p. ± p^ | with the curve r = 1 are possible. But table HI forbids 
two of these. Therefore, only jpg - Pg|, jpg - p^j, and |pg + p^j can occur for disks. 
Yamamoto and Hayashi present experimental evidence in reference 9 to substantiate 
table m for disks. 

A conclusion from this phase of the study is that the maximum number of summed 
and differential nonsynchronous critical speeds is 3 for disks and 6 for pencil- shape ro- 
tors. These results coincide with those found in reference 1 for a firm-foundation model. 

From figure 6 it can be shown that for jig = 1, the curve o> = |p^ + Pg| lies close 
to the curve r = 1 over most of the speed range. Fortunately, table m forbids a non- 
synchronous critical speed for these conditions. Therefore, this potential danger does 
not occur. 


SUMMARY OF RESULTS 

The theoretical critical- speed analysis of rigid rotors on flexible foundations yielded 
these results: 

1. Two frequency modes exist. One is the bouncing mode, which is a function of 
foundation-to-rotor mass and spring- constant ratios only. Another is the conical mode, 
which is a function of moment- of-inertia and spring- constant ratios as well as rotor ro- 
tational speed. 

2. The bouncing- mode natural-frequency result comprises two unique solutions. 

The following results pertain to the conical- mode natural-frequency and critical- 

speed solutions. 

3. Four unique solutions result that are part of low-frequency and high-frequency 
two-branch sets of curves. Forward and backward precession are represented by a 
branch in each set. 

4. Although a firm- foundation model predicts only two of the four critical- speed so- 
lutions predicted by the flexible-foundation model, its two existing solutions agree 

well with their flexible- foundation counterparts. The two omitted critical- speed solu- 
tions can be approximated by using plus and minus the firm- foundation bouncing- mode 
solution (rrj = ±1). 

5. At the foundation-to-rotor spring- constant ratio extremities of zero and infinity, 
the general solution can be approximated by a single two-branch set that represents all 
geometries. This set is identical to the results given by a firm-foundation model. 

6. Because of the existence of this set, only three calculated points are required to 
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cover the entire range of foundation-to- rotor spring- constant ratio for a specified geom- 
etry. 

7. The frequency magnitude of forward precession exceeds that of backward preces- 
sion at all speeds for all rotor shapes. 

8. Individually, increasing rotor rotational speed and rotor polar-to-diametral 
moment- of-inertia ratio each is associated with increased forward-precession and de- 
creased backward-precession frequencies. 

9. All frequencies are inversely proportional to f 1 X /Ml 2 . 

10. Increase in all frequencies is associated with decreased foundation- to-rotor 
moment- of-inertia ratio and increased spring- constant ratio. 

11. No forward-precession high-frequency critical speeds exist for geometries in 
which the rotor polar-to-diametral moment- of-inertia ratio equals or exceeds the 
critical- speed ratio. 

12. In backward precession, disks have lower critical speeds than pencil- shape rotors. 

13. In summed and differential harmonic oscillation, the number of nonsynchronous 
critical speeds varies from 3 for a disk to 6 for pencil- shape rotors. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 12, 1968, 

122-29-02-20-22. 
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